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Abstract 

We present results of applying the Hamiltonian approach to the mass- 
less Schwinger model. A finite basis is constructed using the strong 
coupling expansion to a very high order. Using exact diagonalization, 
the continuum limit can be reliably approached. This allows to repro- 
q_j duce the analytical results for the ground state energy, as well as the 

D vector and scalar mass gaps to an outstanding precision better than 

io- 6 %. 

PACS numbers: ll.15.Ha, 12.38.Gc 

1 Introduction 

The Schwinger model [ij, i.e. quantum electrodynamics in 1+1 dimensions, 
is the simplest gauge theory. Since its formulation in 1962, it has attracted 
much attention. Notwithstanding its apparent simplicity, its physics is sur- 
prisingly rich and in several aspects resembles much more complex theo- 
ries, in particular quantum chromodynamics (QCD). As such, the Schwinger 
model has become the standard toy model for testing lattice techniques^. 

In particular, it was proposed to use lattice Hamiltonian methods for 
investigation of its properties [n-ti^. Using these techniques, several prop- 
erties of the massless and massive Schwinger model were investigated and 
many interesting results were obtained. 

In this paper, we concentrate on the massless case for the 1-flavour model. 
Our main aim is to show that lattice Hamiltonian methods can yield results 
with unprecedented precision - a few orders of magnitude more precise than 
their previous applications. This indicates that their potential has not been 
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fully explored yet. We concentrate on three simple observables - the ground 
state (GS) energy, the scalar mass gap and the vector mass gap. 

In Sec. [2] we introduce the essentials of the Hamiltonian approach and 
the lattice method used. Sec. [3] presents our results. In Sec. HI we sum- 
marize, compare our results with previous estimates and shortly discuss the 
prospects of the method. 

2 Model setup 

We start with the Hamiltonian of the lattice Schwinger model in the Kogut- 
Susskind (staggered) discretization 0, : 

. M 

K = -^Y1 (f{ny 6(n) (t>{n + 1) - 4>\n + l) e -*( n V(»)) + 

n=l 

M 2 M 

+ m £(-l)V(n)0(«) + ^^L 2 (n), (1) 

n=l n=l 

where 4>(n) is a single-component fermion fielcH, defined on each site of an 
M-site lattice with periodic boundary conditions (i.e. on a circle) and obey- 
ing the anticommutation relations {(jy(n), 4>{ m )} = 5 nm , {(f)(n) , 4>(m)} = 0, 
{^(n), (j)^(m)} = 0. m denotes the fermion (quark) mass. The gauge field 
variable 9{n) is defined on the link between sites n and n + 1 and is related 
to the spatial component of the Abelian vector potential by 9{n) = agA(n) 
(we work in the temporal gauge Aq = 0), where g is the gauge coupling and 
a - lattice spacing. The variable L(n) is related to the electric field E(n) 
by the relation L(n) = E(n)/g and to the gauge field by the commutation 
relations: [9(n), L(m)} = iS nm . The possible values of L{n) are quantized: 
L{n)\l) = l\l), where I = 0, ±1, ±2, .... This implies the following action of 
e ±i8(n) on the basig stateg . e ±tf(n)|j) =\l±l). 

For numerical treatment, it is convenient to perform the Jordan- Wigner 
transformation [24| </>(n) = Y[p<n(^ (j3 (p)) a ~ ( n )i where a l (n) are Pauli ma- 
trices (a^ = a 1 ± ia 2 ). This gives: 



M 

Kw = ^Yl {° + (n)e i6 ^o-{n + 1) + a+(n + l)e~ ie ^a- (n)) + 



n=l 



M 2 M 

n=l n=l 



2 The staggered discretization can be viewed as starting with a two-component fermion 
field on each lattice site and ignoring the upper component on odd sites and the lower 
component on even sites. This solves the fermion doubling problem and yields a well 
defined continuum limit in physical observables. 
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Let us now consider the choice of the basis, which is essential for nu- 
merical investigations. The natural choice is the direct product of Ising 
basis, acted upon by Pauli spin operators, and the ladder space of states 
{|/)}, acted upon by the operator L(n) and the rising and lowering opera- 
tors e^ 8 ^. With an M-site lattice, the dimension of the spin part is 2 M , 
while for the gauge part the basis is infinite-dimensional. Hence, the whole 
basis is infinite-dimensional even on a finite lattice. Clearly, for numerical 
computation, some choice of (a finite number of) basis states has to be made. 
One possibility, motivated by the physics of the problem, is to truncate at 
some finite ±/ max . However, such approach gives prohibitively large bases 
even for moderate lattice sizes M, since the region of physically important 
values of I quickly increases as one approaches the continuum limit, adding 
to the exponential increase of the basis size from the spin part. 

Fortunately, a much better truncation procedure exists, using the strong 
coupling expansion (SCE) 

Rewrite the Hamiltonian in a dimension- 
less form: 

w = ^n m = w + xv, (3) 



with: 



05 

M M 



W ° = £2 E (! + (-!)^ 3 (-)) + E L * (»)> (4) 

^ n=l n=l 

M 

V = Y, U + (n)e ie{n) a-(n + 1) + a + (n + l)e" i9(n) a~ (n)) (5) 



71=1 



and 



x = f3 = 1/aV . (6) 

The continuum limit of the model corresponds to x — >• 00, a — >• 0, M — > 00. 
The SCE parameter x is conventionally denoted f3 in lattice gauge theory 
literature. 

Formally, the operator Wo can be treated as an unperturbed part and 
V as a perturbation. It is easy to see that the ground state of Wo, which 
we will denote by |0), is given by: 

|0> = Utlt ••• It) ® |0000- •• 00), (7) 

i.e. "antiferromagnetic" spin state and no gauge field excitations (L(n) = 
for all sites nY Wo|0) = 0. The perturbation operator V flips two neigh- 
bouring spingj and couples them via a gauge field excitation (flux line) 
L(n) = ±1. For example, for a 4-site lattice: 

^|o> = I ;tt4-> ® looo-i) + 1 44tt> ® 1 0010) + 1 tut) ® |o-ioo> + 1 ttU) ® 1 1000). 

(8) 



3 Hence, all generated basis states have zero total magnetization ^ n <r 3 (n) for the spin 
part. 
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Conventionally, a fermionic excitation on an odd site (subscript O; the sites 
are numbered from right to left, n = 1,2,..., M), i.e. | ^o) (I to) means no 
excitation, i.e. spin alignment corresponding to the ground state) is referred 
to as a quark, whereas an excitation on an even site (subscript E), i.e. ||e) 
is an antiquark. Thus, | IeIo) ® 1 0-1) is a meson (quark-antiquark with a 
directed gauge field link connecting them) and, correspondingly | 4oTe) ( X'|10) 
an antimeson (flux line directed in the other direction). 

The above state V\0) = \l s ) is hence a superposition of all possible 1- 
meson ("scalar" - superscript S) states corresponding to one quark - one 
antiquark excitation. We will call such a superposition a class of states. 
A class of states is translationally invariant, however, since we are working 
in the staggered lattice formalism, a meson can only be translated by an 
even number of sites, if an odd number of sites translation is performed, the 
meson is transformed to an antimeson ("charge conjugation" operation). 
One additional symmetry operation that can be performed on each class of 
states is a change of "helicity" H i.e. the clockwise numbering of lattice sites 
to an anti-clockwise (and taking into account the inequivalence of even and 
odd sites), or vice versa. Any of the following operations, or any combination 
of them - translation by two lattice sites T2, "charge conjugation" C and 
"helicity" H - leaves the class of states \l s ) (and any other class) invariant. 

The construction of the dimensionless Hamiltonian W basis is performed 
in the following way. One starts with the ground state |0) and acts with the 
perturbation operator V on it to obtain a 1-meson class of states |1 ). Then, 
one acts again with V, obtaining a 2-meson class of states and a 0-meson 
class of states (i.e. ground state) etc. In general, a fc-meson class of states, 
acted on with V, can produce several (k + 1)- and (k — l)-meson classes of 
st cites * 

V\kf) = J>|(* + l)f ) + - 1)?>, (9) 

3 3 

where \kf) is i-th fc-meson class of "scalar" states and ctj, bj are coefficients 
that have to be calculated for each state. States belonging to a given class 
of states are all states related by the operations T2, C and H. All generated 
classes of states are eigenstates of the unperturbed operator Wo (diagonal 
elements of W), while states with different numbers of mesons are related 
via off-diagonal matrix elements (kf\V\(k + l)j). 

One can also construct the Hamiltonian Wy starting from the lowest 
"vector" (1-meson) class of states (instead of the 0-meson state |0)): 

I 1 "") = (° + (n)e W ^a-(n + 1) - a + {n + \)e~ ie ^ a~ {n)) |0). 

* n 

(10) 

Acting several times with V on \l v ) allows to obtain the "vector" -states 
Hamiltonian matrix Wy in the \kj) basis. 
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Obviously, the construction of the Hamiltonian basis can proceed indef- 
initely, since each site can be connected to a neighbouring one by gauge 
links (flux lines) with arbitrary |-L(n)| = 0, 1,2, . . .. Therefore, we truncate 
at some finite order N of SCE. We will show that a finite N is enough to 
extract the relevant physics. The value of N needed depends on the lattice 
spacing and grows as the continuum limit is approached. A particular role, 
to be discussed below, is played by closed flux lines, i.e. flux loops. When 
the order of SCE N becomes equal to the number of sites M, the following 
class of states appears: 

|0i_ioop) = nut • • • it)®|iin • • • n)+imt • • • it)®i-i-i-i-i • • • (ii) 

with all spins pointing in the same direction as in the ground state (no 
quarks or antiquarks and hence no mesons), but with a flux line connecting 
all sites and hence forming a flux loop (the sign of '1' in the ladder space part 
of a basis state corresponds to either a clockwise or anticlockwise orientation 
of gauge links circulating around the loop). Increasing N, one encounters 
states with an increasing number of flux loops - Eq. ([lip with '1' replaced 
by an arbitrary N\ oop (starting at SCE order TV = N\ oop M). 

To summarize this part, let us shortly comment on the total size of the 
basis and its scaling with the lattice size. In the adopted approach, the size 
of the spin part of the basis grows much slower than exponentially. Zero 
magnetization sector for an M-site lattice consists of [^) states. Group- 
ing states into classes of states by using symmetry operations T2, C and H 
further decreases the spin part basis size. The total basis size is then deter- 
mined by the value of the maximal allowed number of flux loops N\ oop , i.e. 
by the ratio N/M . Note that some given maximal allowed iVi oop does not 
imply in general that / max = N\ oop + 1 , as can be naively expected^. There- 
fore, it is still unavoidable that increasing the lattice size M and keeping 
fixed the ratio N/M, the total basis size grows (approximately) exponen- 
tially. We will give some explicit values of the total basis size in the next 
section. 

In this paper, we are interested in three observables: the ground state 
energy, the scalar mass gap and the vector mass gap. The i-th eigenvalue 
(i = 0, 1, . . .) of the Hamiltonian W will be denoted by Wj. We define the 
following observables (and also include the exact values in the continuum 
limit, for the massless case m = 0) : 



GS energy: 



E = > —» -0.3183098862, (12) 

U 2Mx M^oo 7T V ' 



4 This is most easily understood by looking e.g. at the following state: | t4-ttt444) ® 
1 0-10-1-1-2-1-1), which does not have closed flux loops (i.e. it occurs even if max. ATi oop = 0), 
but still has I = ±2 in the gauge part. 
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x = p 


ag 


lattice sizes M 


max. N 


dimension of Wy 


2500 


0.02 


8,10,12,14 


150 


3347 


10000 


0.01 


8,10,12,14 


200 


4454 


40000 


0.005 


8,10, 12,14 


250 


5561 


250000 


0.002 


8,10,12,14 


350 


7776 


1000000 


0.001 


8,10, 12,14 


600 


13304 


4000000 


0.0005 


8,10,12,14 


1200 


26597 



Table 1: Simulation parameters - inverse coupling x, coupling in lattice units 
ag, lattice sizes M, maximum order of SCE N. We also give the dimension 
of the "vector" Hamiltonian matrix Wy for our largest lattice size M = 14. 

• scalar mass gap: 

M Q UJ-\ — UJo 2 

— = " > -= « 1.1283791668, 13 

g 2^X M^OO yJlT 

• vector mass gap: 

My UjY — Ldn a ^n 1 

— = _ > -= « 0.5641895836, (14) 

g M^oo yJ-K 

where uJq stands for the lowest eigenvalue of the "vector" Hamiltonian 
Wy. 



3 Results 

The parameters of lattices used for this work are presented in Tab. [TJ We 
work very close to the continuum limit and move towards it by changing the 
inverse coupling x = (3 from 2500 to 40000000. The lattices have M = 8 
to 14 sites. The order of SCE (denoted by N) is chosen such that its fur- 
ther increase does not change the results up to machine precision. We also 
illustrate how the size of the Hamiltonian matrix Wy grows with increasing 
N, at fixed M, giving examples for our largest lattice size M = 14. How- 
ever, as we have mentioned above, the size of the Hamiltonian matrix grows 
roughly exponentially when one increases M at a fixed ratio N/M - e.g. 
for N/M = 1200/14, the dimension of Wy is: 1546, 3603, 9615, 26597 for 
M = 8, 10, 12, 14, respectively. 

We first discuss the dependence of the results on the order of SCE N 
and the maximal allowed number of flux loops taken into account in the 

5 Note that typical values of j3 in Monte Carlo simulations for the Schwinger model 
are O(5-10). In the Hamiltonian formalism, one can work much closer to the continuum 
limit. 
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N 

Figure 1: Eigenvalue flow with the order of strong coupling expansion N. 
The vertical scale is split into three regions, for better presentation. M = 8, 
x = (3 = 2500. The difference between the first excited state energy and the 
ground state energy is the scalar mass gap. 

computation N/M. In Fig. [H we illustrate the flow of eigenvalues with N 
for a small lattice with 8 sites and a rather large coupling 1/y/x = 0.02. 
For better visibility, we split the vertical scale into three regions. The most 
distinctive feature of the plot is that the magnitude of i-th eigenvalue is 
saturated for some finite value of N which we denote by iVj. We observe 
that Ni increases as i is increased. Hence, the GS energy is quickly saturated 
at order N\ (only a few flux loops are needed), while the scalar mass gap 
is saturated at a larger value N2. Observables involving higher eigenvalues 
will, of course, require even higher values of N. 

This is further illustrated in Fig. [21 where we show how the three ob- 
servables of interest are saturated with increasing N. We compare different 
lattice sizes and different values of the gauge coupling. It is clear from these 
plots that the GS energy saturates very rapidly, while the convergence is not 
so fast for the mass gaps (it is slightly faster for the vector mass gap). The 
observed convergence to the N — > 00 value is approximately exponential - 
see the lower left plot of Fig. [5] for the scalar mass gap case. This plot shows 
that the saturation order N for some requested precision is highly depen- 
dent on x and M. The contribution of states with an increasing number 
of allowed flux loops becomes more important as one moves towards the 
continuum limit (increasing x). We also observe that at fixed coupling x, 
the saturation order N grows approximately linearly with M, motivating 
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Figure 2: Dependence of scalar mass gap (upper left), vector mass gap 
(upper right) and ground state energy (lower right). We also show (lower 
left) the convergence of the scalar mass gap to the asymptotic (N — > oo) 
value. The dotted lines correspond to M = 8 and solid lines to M = 12. 
Two values of the gauge coupling x = [3 = 2500 (red) and x = (3 = 1000000 
(blue) are used. The exact continuum results (black dotted lines) are also 
shown. 



keeping the ratio N/M, i.e. the maximal allowed number of flux loops, fixed 
when increasing M at fixed x. 

We now discuss the extrapolation of our finite-lattice results to the in- 



finite volume and continuum limits. It was suggested in Ref. [ll| that the 



approach to the M — > oo limit may be performed using the finite size scaling 



method of Fisher and Barber [25J]. However, it was later realized that lat- 
tice data shows rather a power law behaviour close to the continuum limit 
In our data, we clearly find that this is indeed the case. The origin 
of this behaviouin was attributed 18y to the fact that the excitations have 
finite momentum 0(tt/M) and their energies thus obtain a kinetic energy 
correction term 0((7r/M) 2 ). Approaching the continuum limit x — > oo, the 
power law corrections become increasingly important. 

We show our results for the extrapolation to the infinite volume limit 



6 Similar power law behaviour in small-volume QCD was observed in Refs. [2y| and [27 
and models were proposed therein to explain it. 
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Figure 3: The main plot presents the dependence of our lattice data (at 
different couplings x) for the GS energy Eq on 1/M 2 , showing the extrapo- 
lation to the infinite volume limit according to Eq. (|15p . The left inset is the 
close-up of the main plot in the region of M — > oo. The right inset shows 
the extrapolation of infinite volume results at different couplings l/\/x ac- 
cording to Eq. (|16p . Circles at the origin of each plot show the exact infinite 
volume continuum result. Note that the point corresponding to x = 4000000 
is plotted in the right inset, but it is not included in the continuum fit. 



in the main plots of Fig. [3] (for the GS energy Eq), Fig. H] (for the scalar 
mass gap M$/g) and Fig. [5] (for the vector mass gap Myjg). The left 
insets show the neighbourhood of 1/M 2 — > 0. In these plots, the circle at 
1/M 2 = indicates the exact infinite volume continuum result. We fitted 
our data using the following polynomial ansatz (cubic in 1/M 2 ) with fitting 
parameters Fo(ag)\M=oo, F 2 , F4 and Fq : 

F(l/M 2 )\ ag = F (ag)\ M=OD + F 2 ■ ± + F A ■ ^ + F 6 • ^, (15) 

where F(l/M 2 )\ ag stands for either of the three considered quantities: Eq, 
Ms/g or My I g at a fixed coupling ag and for lattice size M. Thus, the 
fitting parameter Fo(ag)\M=oo is the infinite volume result at fixed coupling 
(i.e. at non-zero lattice spacing). We also tried restricting the polynomial to 
a linear or quadratic one in 1 /M 2 . Obviously, the cubic fit has 4 parameters 
and we fit 4 data points - hence the "fits" go exactly through all points. 
However, our precision is such that the inclusion of the quadratic and then 
the cubic term systematically improves our continuum limit value estimate. 
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0.005 0.01 0.015 0.02 0.025 _ 

iNx 



0.01 
2 



1/M 



Figure 4: As Fig. [3l but for the scalar mass gap Ms/g. 




Figure 5: As Fig. [3l but for the vector mass gap My jg. 



Visually, the "fits" are straight lines in 1/M 2 - the quadratic and cubic 
corrections are small, but still relevant from the point of extrapolation to 
the continuum limit. Figs. O H] and [5] show that the infinite volume limit of 
each series of volumes approaches the continuum limit value as x is increased. 
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E M s /g 


M v /g 


exact 


-0.3183098862 1.1283791668 


0.5641895836 




QUADRATIC FITS IN l/M 2 


Fqo 


-0.3183103983 1.1283792781 


0.5641898069 


Foi 


0.009812 -0.013058 


-0.004361 


F02 


0.121528 




error 


1.6-10" 4 % 9.8-10" 6 % 


4.0 • 10" 5 % 




CUBIC FITS IN l/M 2 






-0.3183098827 1.1283791671 


0.5641895845 


Foi 


0.00793666 -0.010565 


-0.003522 


F02 


0.126465 -0.000144 




error 


1.1 • 10" 6 % 2.9 • 10" 8 % 


1.8 • 10~ 7 % 



Table 2: Estimates of fitting coefficients in Eq. (fl~6|) . We quote values for 2 
cases: the infinite volume limit values Fq are extracted either with quadratic 
or cubic fits in l/M 2 . The values of i*bo|cont.,M=oo (denoted by Foo for 
shortness) are our infinite volume, continuum limit estimates and should 
be compared with exact values. The error is K-Foo — exact ) /exact |. The 
quadratic corrections in ag (.F02) are statistically insignificant for the mass 
gaps (and hence not given), except for the case of Ms/g with infinite volume 
taken including a cubic term in l/M 2 . 

The infinite volume values Fo(ag)\M=oo for different couplings ag = 1/V# 
are then plotted against ag on the right insets of Figs. OS] and [5j The fitting 
ansatz for the continuum limit ag — > extrapolation is: 

cont., M=oo 

+ F 01 -(ag)+F 02 -(ag) 2 , (16) 

where the fitting parameters are i r oo|cont.,Af=oo) Fqi, Fq 2 - In this way, the 
extracted values of Fool cont., M=oo are the infinite volume, continuum results 
that can be compared to analytical formulae. The fitting curves for the mass 
gaps, shown in Figs. 0] and El are to a very good approximation straight 
lines in ag, whereas the curve for the GS energy, Fig. [3l shows some mild 
curvature. 

The parameters of these fits are given in Tab. [2j Indeed, the quadratic 
correction Fq 2 ■ (ag) 2 is sizable only for Eq. We observe systematic improve- 
ment of the compatibility of our final continuum limit results i ? oo| C ont.,M=oo 
with the exact result when the term cubic in l/M 2 is introduced in the 
fitting ansatz of Eq. (|15p . It constitutes a very small correction, but at this 
level of precision it improves the infinite volume estimates Fo(ag)\M=<x such 
that the final continuum estimates Foo|cont. M=oo are closer to the exact re- 
sult by 1-2 orders of magnitud^. The best precision is reached for the scalar 

introducing the quadratic term in l/M 2 (in Eq. (I15[) ) is, of course, even more impor- 
tant, improving the continuum limit result by 2-3 orders of magnitude. 
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M s /g M v /g 





result 


error 


result 


error 


exact 


1.1283791668 




0.5641895836 




this work 


1.1283791671 


2.9 • 10"* % 


0.5641895845 


1.8 • 10~ 7 % 


Ref. [11] 


1.120 


0.7% 


0.560 


0.7% 


Ref. [14] 


1.128 


0.03 % 


0.565 


0.1% 


Ref. [18] (I) 


1.25 


11% 


0.56 


0.7% 


Ref. [18] (II) 


1.14 


1% 


0.57 


1% 


Ref. [20] (I) 


1.11 


1.6% 


0.563 


0.2% 


Ref. [20] (II) 


1.1284 


0.002 % 


0.56417 


0.003 % 


Ref. [22] 






0.56419 


7 • 10~ 5 % 



Table 3: Comparison of results for the scalar and vector mass gaps presented 
in this work and in selected literature. 



mass gap, where our estimate is only 2.9 • 10 -8 % off the exact value, i.e. the 
difference is on the tenth decimal place. The results for the vector mass gap 
and GS energy are only slightly worse. 

4 Summary and discussion 

It is clear that the precision of the final continuum limit estimates can be 
arbitrarily increased, by taking larger lattice sizes (e.g. the inclusion of larger 
lattice sizes would allow to introduce a term quartic in 1/M 2 in Eq. (|15p ). 
This would also make it possible to observe sensitivity of the continuum limit 
fit of Eq. (fl~6j) to a term quadratic or cubic in ag (for the mass gaps or GS 
energy, respectively) and perform simulations even closer to the continuum 
limit. 

Obviously, since exact results are known, this is rather pointlesa^]. We 
believe, though, that the results presented in this paper are very interesting. 
Exact results are approached with remarkable rapidity using only moder- 
ate lattice sizes - 8 to 14 lattice sites. The success of the method should 
therefore be attributed to working very close to the continuum limit, at 
inverse coupling x = (3 a few orders of magnitude larger than in typical 
lattice simulations, using e.g. Monte Carlo methods. In earlier applications 
of Hamiltonian techniques, it was also not realized that the continuum limit 
extrapolation can be performed using very high values of x. Using rather 
small values of x, exact results could be approached with a precision of the 

8 Note also that at our current level of precision for the observables, one is approaching 
the problem of machine precision, since typical algorithms for finding eigenvalues can give 
their magnitudes accurate to 0(1CF 10 ) with double precision numbers (several iterations 
of an algorithm with many operations per iteration, accumulating 0(1CP 15 ) errors per 
operation). Hence, increasing precision beyond what we present in this paper would need 
an increased machine precision. 
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order of 0.1-1 % (see Tab. [3]h The impressive improvement of this precision 
with the DMRG method [22] by 2-3 orders of magnitud^l seemed to set 
impassable limits for lattice Hamiltonian methods. However, as we have 
shown in this paper, this precision can be easily surpassed by another 3-4 
orders of magnitude, applying only very simple ideas. 

There is no straightforward way to extend the method presented in this 
paper to the massive case. The reason for this is the following. The mass 
gap Mfij ' g (/i = V, S) at an arbitrary quark mass m/g can be expressed as: 

M/j, (m\ _ M*jp Mn fm 



9 J 9 9 \9 



+ S— - > (17) 



where the first term jg is mass independent and given by the leading 



order term of mass perturbation theory of mass gap expansion 2814301]. i.e. 

or 2/^/7r for the vector and scalar case, respectively. The second 
term SM^/g (m/g), the mass gap shift from massive quarks, depends on 
m/g and has been calculated in mass perturbation theory up to second 
order in m/g. Our data indicate that M^/g can be calculated accurately 
even if the physical system size is very small with respect to the size of the 
meson under consideration (which is the case in our analysis very close to 
the continuum limit, i.e. at very small lattice spacing, and with our rather 
small lattice sizes). However, the shift dM^/g (m/g), incorporating effects 
of quark masses, requires that the physical system size is large with respect 
to meson size. Our small physical system sizes are therefore large enough 
to extract Mj^'/g and finite size effects for this quantity are accurately 



described by Eq. (fT5|) . They are not enough to describe the effects of massive 
quarks, though. This indicates that the nature of finite size effects is different 
for the two terms in Eq. ()17p . Indeed, a naive application of the method 
described in this paper for m/g 7^ yields values for the mass gaps very 
close to the ones in the massless case, i.e. /g is reproduced correctly, 
but bM^/g (m/g) is not. 

A possible solution is then to keep away from the continuum limit (work 
at much smaller x), such that the lattice spacing is coarse enough and the 
physical extent of the lattice is large with respect to the meson size. This 
allows to truncate SCE at some small order N ~ M (with no or very few 
flux loops), which is then enough to saturate the ground state energy and 
mass gaps. This, however, makes it necessary to perform a long extrapo- 
lation to the continuum limit and the amazing precision of the final result 
is lost. Alternatively, to retain this precision, one can simulate close to the 
continuum limit, but using large lattices, such that their physical sizes are 



9 In Tab. [31 the error of the final result for the vector mass gap is 7 ■ 1CT 5 %. However, 
taking into account the estimated error, i.e. Mv/g = 0.56419(4), the actual precision 
claimed by the Authors of Ref. Q is of O(10" 3 ) %. 
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large enough. Our analysis of the massive case will be presented elsewhere. 
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